Nuclear shadowing in Glauber-Gribov theory with (Revolution 



o 

(N 



in 



Nestor Armesto, 1 Alexei B. Kaidalov, 2 
Carlos A. Salgado, 1 and Konrad Tywoniuk 1 

1 Departamento de Ffsica de Partfculas and IGFAE, 

Universidade de Santiago de Compostela, 
15706 Santiago de Compostela, Galicia, Spain 

2 Institute of Theoretical and Experimental Physics, 

117259 Moscow, Russia 

March 16, 2010 



C"| Abstract 
Ph 

We consider deep inelastic scattering off nuclei in the Regge limit within the Glauber-Gribov 
model. Using unitarized parton distribution functions for the proton, we find sizeable shadowing 
effects on the nuclear total and longitudinal structure functions, and F£, in the low- a; limit. 
Extending a fan-diagram analysis for the large-mass region of coherent diffraction off nuclei to 
high Q 2 , we also find significant shadowing effects in this kinematical regime. Finally, we discuss 
shortcomings of our approach and possible extensions of the model to other kinematical regimes. 

l^j. PACS: 12.38.-t, 24.85. +p, 25.30.-c 
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^ 1 Introduction 

O 

The Regge limit of QCD (s — > oo at constant Q 2 ) is of great interest both from the perspective 
of high-energy collider experiments and inasmuch as it poses a challenge for the theoretical under- 
standing of the hadronic wave-function. The rapid growth of cross sections observed at high Q 2 
^ and fairly low x is expected to slow down due to unitarity as x — > 0. Compared to the nucleon 

case, in 7* A collisions these effects are enhanced by the nuclear thickness factor ~ A 1 / 3 . While 
a fully consistent description of nucleus-nucleus collisions at high energies is yet to be found, the 
characteristic unitarity corrections are expected to emerge in the context of deep-inelastic scattering 
(DIS) off nuclei at high energies. 

Nuclear modifications of the total 7* A cross sections are seen to arise for all x > W- 2 0. 
Nuclear parton distribution functions (PDFs) have been extracted at LO and NLO from the avail- 
able data in [ 2—4] . The resulting PDFs differ significantly at low momentum scales due to different 
choices of functional forms and data samples made by the various groups, especially in the gluon 
sector, but gives rather consistent results at higher Q 2 . Yet, the restricted kinematic access of 
present-day experiments give weak constraints on the low-x region interesting for high-energy ex- 
periments and from the point of view of unitarity corrections. 
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The depletion of the nuclear structure function observed at x < 0.1, compared to the incoherent 
superposition of A ^*p cross sections, is usually called shadowing (for a recent review, see [5]). As 
mentioned before, these effects are expected to occur at high energies due to the unitarization of the 
total cross section. In the infinite momentum frame (of the fast nucleus) the depletion is interpreted 
as destructive interference between gluons originating from different nucleons [6]. Similarly, in a 
frame where the target nucleus is at rest, the same effects arise due to multiple scattering of the 
projectile. 

In the latter context, there emerges a critical length scale related to a change of the underlying 
space-time picture of the collision. The coherence length (or life-time) of a given fluctuation of the 
incoming projectile is given by 

l c = \_^AB „ 1 (1) 

Q Q 2mjvx 

in the limit Im^E^AB = W 2 3> Q 2 , where x is the Bjorken variable. At low energies, where lc 
is of the order of the internucleon distance, the projectile undergoes incoherent multiple scattering 
off the target. Remarkably, all higher-order rescatterings cancel and the total "f*A cross section is 
simply given as a superposition of j*p collisions. The critical value is reached when the coherence 
length becomes of the order of the nuclear radius. For lc > Ra, i.e. at x < l/2mjyRA, the 
projectile scatters coherently off all constituents of the nucleus at some given impact parameter. 
Despite the non-local nature of the interactions, the total cross section can be written in the form 
of a multiple scattering series, now including corrections from higher-order rescattering diagrams 
which lead to an overall depletion of the total cross section. 

The analysis of rescattering of the lowest order Fock-states can readily be performed within the 
framework of the dipole model [7j|9]. The dynamics of the scattering is encoded in the dipole-target 
cross section. Whereas the -f*p cross section can be analyzed to great detail due to the wealth of 
available experimental data, the additional impact parameter dependence arising naturally in the 
nuclear case is largely unconstrained. There are several propositions for such an ansatz in the 



literature 10-12 . An additional challenge to this type of models is to properly take into account 



scattering of higher-order Fock-states which dictate the strength of the rescattering and could, 



thus, underestimate nuclear shadowing. This problem was addressed in 13 where the propagation 
of a given massive dipole Fock-state in a large nucleus was modeled and shadowing from qq+Ng 
configurations was included. 

The large probability of rescattering and, thus, a large probability of diffraction dissociation, is 
known to arise from asymmetric dipole pairs and large dipole fluctuations. These involve a strong 
soft contribution already at leading twist in the cross section. A suitable framework, although 
not rigorously established within QCD, to treat these configurations is provided by the reggeon 
calculus [14], where rescattering of the projectile wave function is accounted for by including the 
full mass-spectrum of hadronic intermediate states. This is the so-called Gribov inelastic shadowing 
[15]|16| (for a comprehensive review see e.g. [17] ). 



In a previous paper |18], we have analyzed the total, longitudinal and diffractive ^*p cross 



sections in an unitarized model 19 20 with Q 2 -evolution at leading order (LO) in the strong 



coupling constant. The non-perturbative model [19,20 included unitarity corrections by summing 
all multi-pomeron exchanges, thus providing initial conditions for QCD evolution equations down 
to x ~ 10~ 8 . The aim of the present work is to extend this analysis to DIS off nuclei. This will 
be done in the framework of the Glauber model [2 11 with its extension to the relativistic case due 



to Gribov 15,16 . This formalism relates the inclusive and diffractive j*p cross sections to the 
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corresponding 7* A ones by means of the AGK cutting rules 22 . The original, non-perturbative 
model (l9||20| was thus used to predict nuclear shadowing at low-Q 2 in good agreement with existing 
experimental data |23|. 



Shadowing from Glauber- Gribov theory has also been calculated previously in a similar model 



based on parameterizations of existing data on inclusive and diffractive cross sections [24 - 26 
In particular, using factorization theorems at high Q 2 , shadowing for each parton species was 
calculated. Calculation of gluon shadowing at low x within the same framework was also presented 
in [27]. Nonetheless, the extrapolation of these results beyond the present experimental range, e.g. 
x < 10 -4 , is much less sound since unitarity can be violated for the gluons already at x ~ 10" 4 at 
low momentum scales [28] . 

The model [18] used as basis for our present calculations is free of these shortcomings, and 
can safely be extrapolated down to very low x in the nuclear case too. Thanks to the inclusion of 
scaling violations coming from QCD evolution it can also be used at high Q 2 , a regime relevant 
for future hadronic colliders. While our previous work 18 meticulously included extensions of the 
original non-perturbative formulation to the region of x — > 1 for both the inclusive and diffractive 
cross sections, we will here focus solely on the low-x regime where the summation of intermediate 
large-mass diffractive states gives the leading contribution. Since the present analysis contains only 
contributions from triple-pomeron diagrams, it is therefore valid only for describing large mass 
diffraction - additional ingredients would be required to properly describe low-mass diffraction and 
elastic scattering. Also, our model does not include anti-shadowing, see e.g. 29 30 , which is clearly 
seen in data for 0.1 < x < 0.3. 

The scope of this paper is as follows. In Sec. [2] we will describe the total 7* A cross section as 
a multiple scattering series and re-sum it in the context of two classes of diagrams, the fan and 
eikonal ones. Results for shadowing at low x and relatively high Q 2 are presented and compared 
with experimental data. The comparison to data at smaller Q 2 can be found in |23| . We also 
calculate the nuclear modifications to the longitudinal structure function. In Sec. [3] we calculate 
the cross section of large mass diffraction off nuclei by summing up all fan diagrams. We note that 
Regge factorization, which is seen to hold for diffractive DIS off protons, is violated in the nuclear 
case. We also calculate the ratio of diffractive and total "f*A cross sections as a function of the 
mass of the system. Finally, we present our conclusions in Sec. [4] 
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We proceed by calculating the total 7* A cross section and, thereby, the nuclear structure function, 
F^, assuming that the nucleus is made of independent nucleons. In the coherent regime, despite the 
complex space-time picture of the interaction, the j*A cross section can be expanded in a multiple 
scattering series, such as 



(!) , (2) 



(2) 



Whereas the first rescattering cross section is simply Ao~ f *N, the second and higher-order rescat- 
tering cross sections involve inelastic intermediate partonic states, and can be found using the 
AGK cutting rules 22 . Assuming purely imaginary amplitudes, the second rescattering is simply 



proportional to the inclusive diffractive cross section, given by 



a 



(2) 



/r M max 
d 2 bT 2 A (b) / ™dM 



dM 2 dt 



F 2 , (t r , 



(3) 



t=o 
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Figure 1: A generic diagram in DIS off a nucleus in the Schwimmer model. Each zigzag line is 
considered as a full 7*p amplitude. 



where the nuclear form factor is 

F A (t min ) = j d 2 b,J (bV-tmin) T A (b) , 

and the minimal momentum transfer to the nucleus is t m i n ~ —mj^Xp. The nuclear density, 
PA(b, z), is taken to be a 2-parameter Woods-Saxon normalized to unity, with = J dzpA(b, z). 

The t-integrated diffractive cross section is obtained from the total reduced diffractive cross section 
assuming an exponential i-behaviour with slope Bq, and is given by 



dM 2 dt 



47r 2 a em (B - a'jp In l/ x F ) ^_ ^(3) 



t=o q 2 (q 2 + m 2 ) 



■pF£>(Q 2 ,x f ,0) , (4) 



where the exponential slope Bq is taken from [3l] to be 5.5 GeV 2 . Finally, the double rescattering 
contribution for the deuteron takes the form 



, ; -_ , dt dM 2 

7-oo J Ml. dM 2 dt 



F D (t) (5) 



t=o 



where i*b(i) = e at , a = 40 GeV -2 . The limits on the integral of masses in eq. ^ are = 4m 2 

and M^ ax is found from the condition that Xp ax = O.lQ which guarantees the appearance of a 
measurable rapidity gap in the event. 

The summation of the full rescattering series in eq. ^ is model dependent. At low x we neglect 
effects related to isospin. Here we include two unitarization schemes. Summing all fan diagrams 
we arrive at the Schwimmer model |32| , valid for collisions of a small projectile on a large target 
(see the Appendix for more details). Finally, the total 7* A cross section is given by 

„Sch _ A _ f ,2. T A(b) , . 

a TA - Aa rp Jdb i + {A _ i)TA{b) Q2)/i?2(2 , Q2) • W 



1 Where x P = (M 2 + Q 2 ) / (W 2 + Q 2 
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Figure 2: Nuclear shadowing compared to experimental data from 33 . Solid (black) curves are 



calculated in the Schwimmer model, dashed (orange) ones in the eikonal model (note that the lines 
overlap). 



On the other hand, the summation of multiple rescatterings in the eikonal unitarization gives 

eik_ A f ,2, 1 ~ ex P [-2(A - l)T A (b) gjxp*, Q 2 )/F 2 (x, Q 2 )] 
a rA - Aa rp J 2{A _ i) g2)/i?2(X) Q2) • (7) 

The eikonal resummation gives stronger nuclear shadowing than the Schwimmer model, which can 
be checked by comparing the second non-trivial rescattering term. In eqs. ^ and Q the function 
g{£,,Q 2 ) is defined as 

/(£, Q 2 ) =in j <d ^ B{ X]P ) xjpF^{xjp, f3, Q 2 ) F 2 (t mm ) , (8) 

where B(xjp) = Bq — a'jp In l/xjp. We define the shadowing ratio for the nuclear total cross section 
and, thus, the nuclear structure function, for two nuclei A and B, as Rf 2 {A/B) = Ba^* a/ 1 Aa^* p = 
BF 2 A ( x ,Q2)/AF 2 B (x,Qi)f\ 

Unfortunately, there is little available data for the kinematics where our model is valid, i.e. 
x < 1CT 2 and Q 2 > 2 GeV 2 . Utilizing the results on X]pF$(xp, /3, Q 2 ) and F 2 (x, Q 2 ) found in 
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we compare shadowing calculated in the Glauber-Gribov theory (the curves for the Schwimmer and 
eikonal unitarizations overlap) to experimental lepton-nucleus data [33] in Fig. |2j Although the 
model calculations are connected by solid lines, we have made calculations at the corresponding 
(Q 2 ) of each data point. Furthermore, we check the Q 2 -dependence of shadowing in Fig. [3] compared 
to experimental data from [34]. Shadowing for a lead (Pb) nucleus at several momentum scales is 
shown in Fig. [4] (left), and the ^-dependence is presented the right part of Fig. [4j 

A comparison to several models of nuclear modifications is presented in Fig. [5] Our model is 
in good agreement with the recent calculation within the dipole model calculation of |13| , where 
rescatterings of the full qq+Ng fluctuation is taken into account, see left plots in Fig.[5j The higher 
Fock-states of the dipole correspond to the summation of triple-pomeron diagrams in our approach. 

Furthermore, we obtain reasonable consistency at low x between our result and the calculations 
presented in [26] , right part of Fig. [5j Different from our work, in [26] the input partonic densities 



2 We automatically define Rf 2 (A) = F^(x, Q 2 ) J ' AF?.(x, Q 2 ), where F2 is the structure function of a proton. 
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Figure 3: Nuclear shadowing compared to experimental data from [34] . Solid (black) curves are 
calculated in the Schwimmer model, dashed (orange) ones in the eikonal model (note that the lines 
overlap). 




Figure 4: Q 2 -dependence of nuclear shadowing for lead (Pb) (left) and shadowing for several lighter 
nuclei at Q 2 = 3.5 GeV 2 (right). Upper (black) curves are calculated in the Schwimmer model, 
lower (orange) ones are calculated in the eikonal model. 
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Figure 5: Results for ^(Pb)/^ compared to the dipole calculation in 13 (left), the EPS09 
NLO parameterization [4] (right) and the FGS model [26] (both). Solid upper (black) curves are 
calculated in the Schwimmer model, lower (orange) ones are calculated in the eikonal model. 



of the proton and pomeron were taken from perturbative fits of experimental data at NLO. Other 
differences stem mainly from the treatment of the large- x region, e.g. the parameter xp max is 



put to 0.03 in 26 and strong anti-shadowing is included in the parameterization of the gluons 
stemming from the momentum sum rule of the nucleus. The latter corresponds to secondary reggeon 
exchanges [301 which are neglected in our model. The agreement at high-Q 2 is also reasonable. 
Finally, we compare our calculations to the recent parameterization of shadowing extracted from 
a NLO DGLAP analysis of nuclear data [4] in the right part of Fig. [5j The errors of the nuclear 
parameterization are also indicated. We find large differences between the approaches at low 
momentum scales at x < 10~ 2 , which are solely due to the choice of initial conditions. The model 
in [5] predicts a significantly smaller nuclear shadowing at x < 10 -4 . The discrepancy persists at 
higher Q 2 . 

It is worth pointing out the different treatment of Q 2 evolution between the different models. 
Our treatment of the Glauber-Gribov rescattering series is valid for both soft and hard processes 
at low-x and involves effectively all twist contributions. Nevertheless, it was demonstrated for the 
non-perturbative regime in [231 that the leading- twist contribution is by far the strongest. Note the 



difference with 26 , where the nuclear PDFs at some initial scale were used as the input to DGLAP 
evolution. The same approach was applied in (4j, where initial nuclear PDFs were fitted. Finally, 
the dipole model calculation [13] also account for some scaling violations. The Q 2 evolution of the 
nuclear structure function is clearly logarithmic in all the discussed models. 

We continue our study by analyzing the longitudinal part of the total cross section. In the 
naive quark-parton model, the Callan-Gross relation gives that Fl = 0. This is not the case when 
we include radiation of quarks and gluons. The longitudinal structure function can readily be 



calculated using the Altarelli-Martinelli formula 35 



relating Fl quantities at LO on the right-hand side of the equation, where Nj is the number 
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Figure 6: Shadowing of the nuclear F 2 (dotted) and Fl (solid) structure functions together with the 
gluon nuclear distribution (dashed curves). Upper (black) curves are calculated in the Schwimmer 
model, lower (orange) ones are calculated in the eikonal model. 



of active quark flavours. Due to large magnitude of the gluon PDF compared to quarks in the 
proton, this quantity is dominated by the gluon component. The corresponding relation for the 
nuclear case is analogously found, so that we can define the nuclear suppression factor, Rp L = 
F^(x,Q 2 )/ AFl(x,Q 2 ), of the longitudinal structure function as 

£ f (f) 2 (Efee? (i - i) yR g (y,Q 2 )g(y,Q 2 ) + lR F2 (y,Q 2 )F 2 ( y ,Q 2 )) 

Rp L = 2 ) (10) 

£ f (§) e * i 1 ~ f ) w(v> ^ + f ^ Q 2 )) 

where we calculate with Nf = 4. In the equation above, R g denotes gluon shadowing and is 
calculated analogously to Rp 2 by substituting f(£,,Q 2 ) by 

/(£, Q 2 ) y 47r f ^ x P f F/p (xp) Pgptf, Q 2 ) F 2 A {t min ) , (11) 

where fp/ p (xp) is the i-integrated pomeron flux factor and gp(fl, Q 2 ) is the pomeron gluon distri- 
bution, and F 2 (x, Q 2 ) — > xg(x, Q 2 ) in eq. Both of these quantities are taken from fl8] . 

Results for Rp L (x,Q 2 ) together with R g (x,Q 2 ) are shown in Fig. [6j Due to the small gluon 
shadowing, nuclear shadowing of the longitudinal structure function is almost the same as for the 
total cross section. The smallness of R g , on the other hand, is caused by the rapid small- x growth 



of the gluon PDF in the proton at LO [18] . Note that studies at NLO 24,25.27 infer a more 
pronounced gluon shadowing, resulting in a significantly smaller Rp L than shown in Fig. |6j see 
e.g. [25]. 



3 Diffraction off nuclei in the Schwimmer model 

Diffraction in DIS provides valuable complementary information to the total cross section as it 
probes the parton distributions at large impact parameters and involve a large soft component. 



8 



Figure 7: The single diffractive cross section in the Schwimmer model. Note that the double zig-zag 
lines at the bottom sums all fan diagram contributions. The top incoming line corresponds to a 
virtual photon, while the lower one to a target nucleon/nucleus. 



Nuclear DIS possess a rich structure and have recently triggered significant activity [9,25, 36-43 1 . 
Generally speaking, one naturally expects a large diffractive cross section in the nuclear case due to 
the large probability of rescattering. Several groups predict a strong enhancement of the fraction 



of diffractive events, up to 30-40% 19 , 40 , 43] (compared to 15% in j*p at HERA), in the nuclear 



case. Note that most of this enhancement derive from the elastic component, while scattering of 



higher-order Fock-states are found to be suppressed 25 41.431. This is in line with the standard 



Good- Walker picture of diffraction, and stems from the fact that as the target gets closer to the 
black disc limit for central and semi-central impact parameters at high energies, only the periphery 
of the nucleus contributes to inelastic diffraction. 

In this work, we focus uniquely on the diffractive processes where a large mass is being produced. 
In this limit, the summation of fan diagrams will dominate the cross section. Although the mass 
spectrum is steeply falling, ~ l/M 2 , the experimental reach of planned lepton-hadron colliders, 



such as e.g. eRHIC 44 and LHeC [45], will be sensitive to these effects. Nuclear diffraction 
at /3 — > 1 is large, and an extension of the model to this region, where rescatterings are very 



important, would be e.g. to include vector-meson dominance, in the spirit of 146,47 . Besides, we 



do not include the possibility of intermediate excitations of target protons, e.g. by introducing a 
two-channel model as in |39|. Additionally, we consider totally coherent diffraction, i.e. assuming 
that the nucleus stays completely intact after the interaction with the projectile. In this case, the 
momentum transfer to the nucleus must be very small, \t\ ~ l/R\- This situation seems more 
easily tractable from the theoretical point of view, whereas the calculation of incoherent diffraction 
involves calculation of survival factors with process dependent absorptive cross sections [39] which 
are a priori unknown. One the other hand, the full i-dependence of diffraction off nuclei would 
have to be modeled including break-up mechanisms and multi-fragmentation. Measuring the intact 
nucleus at small t at electron- ion colliders will be challenging [44]. Efforts to extend the present 
analysis of nuclear diffraction along these lines are under way. 

After resumming all fan diagram exchanges, see Fig. [7] (the derivation can be found in the 
Appendix), the single-diffr active 7* A cross section at a given center-of-mass energy, W 2 (where 
x ~ Q 2 /W 2 + Q 2 ), and scale, Q 2 , as a function of the diffractively produced mass, M 2 , is given 
by |48l|49] 



da? rA 



drP 

7* N 



dM 2 dM 2 



d 2 b- 



4TrB A 2 Ti(b) 



1 + AT A (b) {2f(x™*, Q 2 ) - f(x/f3, Q 2 )) / F 2 (x, Q 2 



(12) 



where /(£, Q ) is defined in eq. d8]) and we have assumed an exponential t-slope of the diffractive 
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Figure 8: Nuclear diffractive ratio in the Schwimmer model with CFSK PDFs. 
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Figure 9: Nuclear diffractive-to-inclusive ratio in the Schwimmer model with CFSK PDFs. 



7*p cross section with slope parameter Bq. Clearly, eq. (12) does not factorize the dependence of 
kinematical variables related to the emission of a (generalized) pomeron from the target (typically 

and t) from those related to its subsequent scattering with the 7*-projectile (typically f3 = 
Q 2 1 ' M 2 + Q 2 and Q 2 )- This indicates that Regge factorization [50] , which is roughly valid in a 
large kinematical range in diffractive DIS off protons |18| , is more strikingly violated for nuclei due 
to the stronger rescattering in the latter case. 

In fact, a peculiar feature can be deduced from eq. (12), or simply by looking at Fig. [7J namely 
that shadowing effects tend to vanish with increasing of the diffractively produced mass. This 
is because, as long as we demand that the nucleus stays intact, the phase space available for 
rescattering is reduced with increasing M 2 . 
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We define the nuclear diffractive ratio in the standard form as 



R v (s,M 2 ,Q 2 



dM 2 



drP 
dM 2 



(13) 



for coherent diffraction^] We plot the suppression factor for Pb as a function of the diffractively 
produced mass for two different momentum transfers in Fig. [8] for three different scenarios: (i) 
"eRHIC low" - assuming a beam energy of Ea = 100 GeV/c per nucleon for the nucleus side and 
an electron beam energy of E e = 10 GeV/c; (ii) "eRHIC high" - the same nucleus beam energy 
as in the previous case and E e = 30 GeV/c; and (iii) "LHeC" - with Ea = 2750 GeV/c per 
nucleon and E e = 100 GeV/c. We stress that this is the suppression solely due to coherent multiple 
rescattering 25,41 . 

In Fig. [9 we plot the ratio a v j a tot for protons and three nuclei ( 12 C, 40 Ca and 208 Pb) for a 
fixed value of /3 and momentum transfer Q 2 as a function of Bjorken x. For such values of /?, where 
we safely can assume the dominance of fan diagrams (M 2 > 135 GeV 2 ), we observe the suspected 
suppression compared to the proton case. It is also worth noticing the different behavior with x of 
the ratio a v j a tot for protons and nuclei at low momentum scales Q 2 . 

Note that we systematically find a nuclear suppression of diffraction, to be contrasted with the 
enhancement found for smaller masses in [9p0l|43]. 



4 Conclusions 

Shadowing of the nuclear structure function was calculated in the Glauber-Gribov model down 
to very low- a? and high-Q 2 utilizing the results on inclusive and diffractive cross sections in j*p 
scattering. The agreement with data is good, and comparison with predictions from other models at 
low-x seems to be fairly consistent, both regarding the magnitude and Q 2 -dependence, except that 
nuclear shadowing of is found to be significantly stronger at low x than in [4] . Strong shadowing 
effects are expected at future high-energy electron-nucleus colliders and heavy-ion collisions at 
LHC 1 17] . On the other hand, weak shadowing of the gluons, owing to the large gluon density of the 



proton at low x and Q 2 found in 18 , lead to significantly smaller nuclear effects on F^, in contrast 
to [25] . Calculating high- mass coherent diffraction off nuclei from fan diagram resummation, we 
found a large suppression solely due to coherent scattering. Also, Regge factorization is strongly 
broken in the nuclear case. 

The program of extending the above j*A analysis to the regime of large x, and low- mass 
diffraction, poses several difficulties that are not yet unambiguously resolved in the literature. 
In the inclusive case, vector- meson dominance and, more seriously, reggeon exchanges with the 
target are not included. While the former can be straightforwardly incorporated (see e.g. |26|), the 
inclusion of the latter is related to anti-shadowing (the second rescattering term is positive) and 
would involve a much more careful treatment of the multiple scattering series. For diffraction off 
nuclei, the most serious challenge is posed by the treatment of low-energy nuclear effects related to 
the break-up of the nucleus - a necessary ingredient for a proper comparison to an experimental 
situation. These lacks can only be amended through a full-fledged Monte-Carlo procedure. An 



estimate of fully-incoherent diffractive scattering was presented in 39 , 43 



3 Note that in the absence of nuclear effects (and assuming a Gaussian density profile for the nucleus) Rt> = 
2B /R 2 A ~ A" 2/3 which is the case when R 2 , — 2Bq. Hence, actually one should compare the quantity A 2 ' 3 Rt> to 
unity. 
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A Derivation of the nuclear total and diffractive cross section in 
the Schwimmer model 



The cross section for diffractive production of a system with large mass in hadronic reactions is 
found by analyzing the contribution from triple-pomeron diagrams. In the situation of asymmetric 
collisions at high energies, e.g. DIS off nucleons and nuclei, diagrams with a single pomeron emitted 
from the projectile with its subsequent branchings give the largest contribution, i.e. that enhanced 
by the product of 3iP-couplings and A 1 / 3 . The resummation of all triple-pomeron branchings, the 
so-called fan diagrams, was first carried out in [32[| and more recently studied in 48,49], where 



expressions for total, inclusive and single-diffractive cross sections can be found together with 
probabilities for large rapidity gap survival. The fan resummation with n — > m pomeron vertices 
was presented in 51 . Here we recall some key results. 

First 
mation. 



, let us define a quantity which will prove helpful in describing the fan-diagram resum- 
Consider scattering off a proton by means of a single triple-pomeron diagram. Assuming 
some maximal allowed mass, or a minimal rapidity gap, yj\j = Y — y = In (M 2 /so), the expression 
is given by 

Y-vm 



°yp(ri*N -»■ XN) 



dyd bd b' g ]PN (b)g J *]p(b 



b')G 3IP e A ye 2 ^ Y - 



■v) 



(14) 



g~i*ipgiPN{b = 0)e AY ■ K N (Y,y M ) 



where we have assumed Gaussian impact parameter profiles. The first part in the expression above 
is the expression for a simple one-pomeron exchange while 

QpnGzp 1 



,A(Y-Vm) 



1 



(15) 



7rAjv A 

where G32P is the triple-pomeron vertex. Due to the smallness of this parameter, found in phe- 
nomenological studies, kn is usually considered negligible. Yet, in the nuclear case it is enhanced 
by the mass number, KA(Y,yM,b) = ATA(b)KN(Y, Vm) (instead of the usual Gaussian parameteri- 
zation, we will use the Woods-Saxon profile for the nuclear case) and is therefore sizeable. On the 
other hand 



M 2 

M 2 dM dM'i dt 



dafjp 



t=o 



where Mi 



4m 2 . In the following we will therefore make the following extension 



KA(Y,y M ,b) oc 



AT A {b) 
<Jip(Y) 



M 2 



dM 



12 



da^p 



Ml 



dM' 2 dt 



(16) 



(17) 



t=o 
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Generalizing this expression to multi-pomeron exchanges, including proper normalization and in- 
cluding coherence effect via the nuclear form factor FA(t m in) [17|, we finally obtain 



K A (Y,y M ,b) 



AitAT A {b) 
(XT 



a 



tot 



M 2 

M 2 dM dM' 2 dt 



da v 



(18) 



t=o 



Sometimes KA(Y,yM,b) is called an effective rescattering cross section. Note that the total cross 
section in the denominator depends on the full rapidity distance, i.e. the center-of-mass collision 
energy. 

Then, the total cross section at a given impact parameter from the summation of all fan 
diagrams is given by [32] 



<r** ch A(Y,b) 



e 9iPr9iPA(b) 
KA (Y,Y,b) + l ' 



(19) 



where the integration limit in ha goes up to a maximal mass, which is usually chosen to be 
constrained by Xp ax = 0.1. The total cross section in the Schwimmer model is given in terms of 
DIS variables in eqs. Q and The total inelastic cross section [48, 49 



is 



inel, Sch / ^r ,\ 
° 7 *A ( y > b ) 



e gprgpAjb) 

2K A (Y,Y,b) + l ' 



(20) 



and, finally, the cross section of single diffraction of the incoming hadron into a hadronic system 
with mass, y = In (M 2 /sq), larger than yM, y > yu, is given by 



(7 



V,Sch 
-y*A 



(Y,y M ,b) 



giP^diPAib) (ka(Y, Y, b) - ka(% vm, b)) 
[2k a (Y, Y, b) - k a (Y, y M ,b) + 1] [k a (Y, Y, b) + 1] ' 



(21) 



Then, taking the derivative, we find the cross section for single diffractive production of a system 
with mass M 2 as 



, V,Sch 
d(T rA 

dyM 



which can be written as 

, V,Sch 
d(7 TA 

dM 2 



d 2 b _ g 2 P A(b)giPTGyp^ {Y+yM) 



[2 KA (Y,Y,b) - K A (Y,y M ,b) + 1Y 



d 2 b- 



AT A {b)ajp(Y) 



[2K A (Y,Y,b) - K A (Y,y M ,b) + l} 2 

^A 2 T\{b) 



dyn 



K A {Y,y M ,b) 



dM 2 dt 



t=o 



d 2 b- 



[2 KA (Y,Y,b) - K A (Y iyM ,b) + IY 



(22) 
(23) 

(24) 



Going to DIS variables, we define the Bjorken-x related to the collision energy, or the full ra- 
pidity distance, Y — > x = Q 2 / (s + Q 2 ) and the mass of the produced system to yu — > xp = 
(M 2 + Q 2 ) J (s + Q 2 ) = x/fi. The final formula is given by eq. (fl2). 

The fan diagram summation presented above was derived in the "soft" limit, i.e. assuming 
simple pomeron exchanges. The generalization to full 7*p cross sections in eq. (18) is an ansatz 
that is known to be valid for the second-order rescattering contribution to the total cross section, 
see eq. This prescription respects QCD scaling violations, and thus we extend it to the full 
analysis of hard reactions off nuclei. 
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